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ABSTRACT 

A study  of  turbulent  diffusion  in  the  surface  boundary  layer 
is  conducted  to  determine  the  effects  of  the  wind  speed  on  the  dis- 
tribution of  sulfur  dioxide  (SO2)  with  time.  With  the  use  of  ■'uto- 
correlograms  and  the  spectra  of  wind  speed  and  SO2.  similarities  be- 
tween these  scalar  quantities  under  varying  stability  conditions  are 
investigated.  Pollution  monitoring  sites  in  the  Salt  Lake  Valley 
vicinity  and  St.  Louis  are  used. 

Under  stable  conditions,  the  auto-correlograms  and  the  spectra 
of  wind  speed  and  SO2  are  similar  in  many  respects.  The  auto-correlo- 
grams show  very  similar  coefficients  up  to  a time  lag  of  13  hours,  and 
the  spectra  shows  similar  peaks  near  periods  of  24,  12  and  8 nours. 

The  wind  speed  peaks  are  more  prominent,  which  suggests  the  spectrum 
of  SO2  is  affected  by  the  spectrum  of  the  wind  speed  for  periods  equal 
to  or  less  than  24  hours.  Forecasting  the  SO2  spectrum  from  the  wind 
speed  spectrum  is  possible  for  the  periods  near  24  and  12  hours. 

The  similarity  between  the  wind  speed  and  SO2  auto-correlograms 
and  spectra  decreases  as  stability  decreases.  In  the  stable  periods, 
the  effect  of  the  large  scale  motions  is  suppressed,  and  more  energy 
is  found  in  the  higher  frequencies  than  in  neutral  or  unstable  periods. 
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Chapter  1 


INTRODUCTION 


In  the  study  of  boundary  layer  diffusion  with  application  to 
air  pollution  concentration  forecasting,  many  forecast  models  have 
been  developed.  The  most  Important  variables  In  most  of  the  models 
are  the  mean  wind  at  varying  levels  through  the  boundary  layer  and 
diffusion  parameters  such  as  turbulence  intensity  and  thermal  strati- 
fication. 

It  is  known  that  the  power  spectrum  of  wind  which  is  a vector 
quantity  is  not  necessarily  similar  to  that  of  pollution  concentra- 
tion which  is  a scalar  quantity.  However,  since  diffusion  is  a con- 
sequence of  turbulent  motion,  the  power  spectrum  of  turbulent  wind 
speed  which  is  a scalar  quantity  may  bear  some  similarity  to  that  of 
pollution  concentration.  Should  such  a similarity  be  established, 
portions  of  the  pollution  concentration  spectrum  can  be  predicted  by 
the  spectrum  of  wind  speed.  The  purpose  of  this  study  is  to  test 
such  an  hypothesis. 

The  autocorrelation  coefficients  and  the  spectra  for  wind 
speed  and  SO2  were  calculated  for  air  pollution  monitoring  sites  in 
the  Salt  Lake  City  vicinity  and  from  urban  and  suburban  sites  in  the 
St.  Louis  area.  Salt  Lake  City  and  St.  Louis  were  chosen  for  their 
different  physical  environments  and  amount  of  data. 
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This  paper  consists  of  five  major  sections.  The  first 


section  deals  with  statistical  theory  which  »s  applicable  to  calcu- 
lating diffusion  parameters.  The  second  section  deals  w1tl\  the 


methods  used  In  this  paper  for  calculating  the  autocorrelation  coef- 
ficients and  spectra.  The  third  section  describes  the  topography  of  I 


the  Salt  Lake  vicinity'  and  St.  Louis,  and  the  type  and  location  of 
monitoring  sites.  The  fourth  section  contains  a discussion  of  how 


the  stability  categories  were  determined  and  the  quantity-  of  data 
available.  In  section  five,  the  results  and  conclusions  of  this 
study  are  presented.  Finally,  a quick  suniiiarv  of  the  conclusions  Is 
In  the  last  section. 


1 


Chapter  2 


r 


APPLICABLE  STATISTICAL  THEORY 


In  the  study  of  atmospheric  diffusion  two  basic  theories 
exist:  the  K or  gradient  transport  theory,  and  the  statistical 
theory.  The  K-theory  is  derived  in  Eulerian  space.  The  values  of 
eddy  diffusivi ties  derived  from  K-theory  appear  in  the  Gaussian  con- 
tinuous point  source  equation  (1)  and  are  a function  of  space.  The 
statistical  theory  is  derived  in  Lagrangian  space  and  the  correspond- 
ing Gaussian  continuous  point  source  equation  is  (2)  in  which  c^(t) 
and  o,(t)  are  a function  of  time.  Other  models  also  employ  the  use 
of  Lagrangian  diffusion  parameters  (c). 


r ? 


C^(x,y,r) 


C ( X ,y , 2 ) 


AnUXyyK^^ 


A. 


exp 


^20y^(t)  2o/(t) 


(1) 

(?) 


<r(x,y,z)  is  the  average  concentration  at  the  desired  point,  u is  the 

mean  wind  velocity,  and  Q is  the  rate  of  emission  of  effluent. 

2 ’’ 

The  values  of  and  K^,  are  related  to  (t)  and  o,‘’(t). 

2 ^ 

As  t-*«'  0 (t)— *2K  t and  c,‘'(t)--*’2K,,t.  Thus,  K-theorv  is  a special 
j yy 

case  of  the  statistical  theory.  But  in  mesoscale  diffusion,  such  as 
the  Salt  Lake  Valley  and  St.  Louis  metropolitan  area,  t is  of  the 
order  of  hours.  Therefore,  statistical  theory  is  more  accurate  for 


L 


these  arens. 

Equation  (2)  is  used  when  calculating  pollution  concentration 
downwind  of  a ooint  source  which  is  emitting  effluent  at  a steady  con- 
tinuous rate  over  flat  terrain.  The  most  important  terms  in  the 
equation,  or  any  pollution  model  equations,  are  the  u and  the  diffu- 
sion oarameters.  u must  be  steady  and  uni-directional  and  o^lt)  and 
c',(t)  are  a function  of  atmospheric  turbulence,  space  and  time.  This 
work  is  aimed  at  a better  understanding  of  the  diffusion  parameters 
(Oy‘'(t)  and  «tid  how  they  vary  with  atmospheric  stability. 

In  order  to  apply  the  Gaussian  equation  to  predict  concentra- 

9 2 

tion,  0 ‘‘(t)  and  o (t)  must  be  known.  G.  I.  Taylor  (1921)  found  the 

y * 

following  relationship: 


0 ‘(t)  - 2/1  <v,'(t)v,'(t')>  dt'dn 
0 0 


where  is  the  Lagranqian  turbulent  velocity,  n is  a dunr\v  variable 
and  vv' (t)v' (f )''  is  the  autocorrelation  function,  which  may  be  cal- 
culated as  an  ensemble  average  of  many  data  samples  called  realira- 
tions. 


If  V|^  is  plotted  as  the  ordinant  and  t as  the  abscissa  for 


all  realizations  then 


wher^  N is  the  number  of  realizations  or  trails.  If  the  turbulence 
is  homogeneous  and  stationary,  (4)  can  be  written  as  <Vj^(t)V|^ (t+T  )> 
where  t is  the  time  laq  defined  as 


5 


T = t'  - t (5) 

Now  (4)  is  a function  of  t not  t.  Tho  autocorrelation  coefficient 
may  then  be  defined  as 


<V||^(t)V|^(t+T)>  V|J^"(t)V|^(t+T) 

''l  <v-^(t)>  vp(0 


(6) 


9 

where  is  used  to  normalize  the  equation.  Equation  (3)  can  now 

be  written  as 


R (T)dTdri 

''l 


(7) 


From  integration  by  parts  of  (7),  Kampe  de  Feriet  (1939)  de- 
rived the  following  equation: 


c/(t) 


2^2 


f (t-T)R„  (T)dT 


(8) 


2 

In  order  to  calculate  (t),  the  integral  in  (8)  must  be  known.  The 
integral  can  be  approximated  by 

OP 

/ R (T)dn  = (9) 

0 L 

where  is  the  Lagrangian  Integral  time  scale.  Equation  (8)  for  long 
diffusion  time  may  be  expressed  as 


(t)  5 2v 


(10) 


Evaluations  of  (i)  and  were  done  for  the  Magna  station 


in  the  Salt  Lake  Valley  and  are  presented  in  Chapter  6,  Results  and 
Conclusions. 

With  the  advent  of  the  Fast  Fourier  Transform  algorithm,  a 
more  useful  and  informative  method  of  determining  the  diffusion  param- 
eters was  found.  This  method  involves  the  use  of  the  power  spectrum, 

S (a')  of  the  Eulerian  wind  speed  v^.  Equation  (11)  shows  the  rela- 

2 

tionship  of  (w)  to  R (t)  and  c (t). 

a ''l  ^ 

O' 

S (a')  “ ^ / R,,  (t)  cos  a'Tdt  (11) 

''L  0 ''L 

With  much  manipulation, 

c./(t)  . d.  (12) 


Equation  (12)  shows  the  diffusion  parameters  can  be  found  bv  knowing 
only  the  spectrum  and  the  diffusion  time  of  interest. 

The  data  that  is  used  in  this  study  was  gathered  at  fixed 
stations  which  is  the  Eulerian  system.  However,  (10)  and  (12)  are  in 
the  Lagrangian  system.  The  Equations  (10)  and  (12)  must  now  be  con- 
verted into  the  Eulerian  system. 

Pasquill  (1961)  found  that  for  the  microscale,  the  Lagrangian 
and  Eulerian  time  scales  are  related  by 

B = s 20  (13) 

'e 

For  large  diffusion  time  in  the  Eulerian  system  (10)  may  be  expressed 


1 


i. 


For  the  wind  speed  and  SO2  we  want  to  determine  where  the 
energy  is  located  and  what  is  the  rate  of  dissipation  in  the  higher 
frequencies.  This  information  will  give  a more  concise  understanding 
of  the  S (w)  in  (15)  and  help  in  forecasting  the  SO,  spectra  from  the 

''e 

wind  speed  spectra  which  will  ultimately  lead  to  improved  forecasting 
of  SO2  concentrations  at  different  locations. 

2 2 

With  equations  (14)  and  (15),  good  values  of  (t)  or  0^  (t) 
should  be  obtained  providing  homogeneous  and  stationary  turbulence 
exist.  These  conditions  do  not  exist  in  the  large  scale  in  either 
the  Salt  Lake  Valley  or  in  St.  Louis,  so  many  realizations  must  be 
taken  under  differing  stability  conditions  to  evaluate  how  the  auto- 
correlation coefficients  and  the  spectra  behave.  With  many  relations, 
quasi -stationary  turbulence  can  be  assumed  for  a given  stability  con- 
dition and  the  statistical  theory  can  be  applied  to  forecasting  pollu- 
tion concentrations. 


CHAPTER  3 


METHODS  OF  CALCULATION 


R and  R. 


The  most  important  variable  in  Taylor's  diffusion  equation  is 
the  autocorrelation  function  R or  R and  how  it  varies  with  chang- 

'^E  ''e 

ing  time  lag.  The  classical  method  for  calculating  R is 

''e 


Jl  u'.(t)u:(t+T) 
, M ~ 

R j,  “r<‘' 


where  M is  the  length  of  the  interval  frjm  which  R^  is  calculated. 
Box  and  Jenkins  (1970)  suggest  the  following  method: 


^ u^(t)uUt+T) 

^ J J 

TN ■ \ H ; 

I I U^^(t)  V Ul"{t+T) 

J J 


where  N=n-T  and  n is  the  nunibe*"  of  observations.  Fioure  1 shows  R 

''e 

calculated  for  the  wind  speed  at  Magna  under  neutral  conditions.  Note 
the  values  of  R^  for  the  classical  method  are  slightly  lower  than 
that  for  the  Box  and  Jenkins  Method,  however,  the  shape  of  the  curve 


is  similar. 

Two  types  of  averages  were  used  in  this  study,  the  24  hour 


moving  average  and  eliminating  the  linear  trend. 

The  24  hour  moving  average  was  calculated  using 

“o  ' “j-12  " 


10 


(18) 


where  and  Uj.12  values  plus  and  minus  12  hours  from  the 

u.  in  question.  A il.  is  calculated  for  every  point  in  the  realization 

J V 

and  u.  is  found  by 

uj  . Uj  - S.  (19) 

The  linear  trend  was  the  other  method  used  to  calculate  i7.. 

w' 

Any  steady  increase  or  decrease  in  concentration  or  wind  speed  which 
prevailed  throughout  the  entire  realization  is  eliminated.  Equation 
(19)  is  then  used  to  calculate  uL 

J 

The  Spectrum 

J.  N.  Rayner  (1971)  outlines  a direct  method  to  calculate  the 
spectrum.  I have  used  this  method  and  modified  it  slightly  to  accomo- 
date spectral  density  and  ensemble  averaging  calculations.  See 
Appendix  F for  computer  programs. 

The  first  step  of  the  direct  method  is  to  remove  the  linear 
trend.  Since  periods  as  short  as  three  days  in  the  Salt  Lake  data 
and  one  day  in  the  St.  Louis  data  are  used,  any  wave  longer  than  three 
days  or  one  day,  respectively,  will  he  affected  by  this  trend  removal 
and  might  not  produce  meaningful  results. 

Once  the  trend  was  removed,  a cosine  bell  window  function  was 
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multiplied  to  the  first  and  the  last  ten  percent  of  the  data  in  the 
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following  manner.  For  the  first  ten  percent  of  the  data,  the  window 
function  (20)  is  used. 


G (20) 

where  h[j]  is  the  window  function,  D is  the  length  of  the  time  series, 
G = D/10,  and  j is  the  data  hour.  For  the  last  ten  percent,  the  de- 
scending cosine  bell  window  function  is 

0 - G < j < D - 1 (21) 

The  third  step  is  to  add  zeros  to  both  ends  of  each  time 
series.  This  is  done  for  two  reasons;  one,  to  correct  for  the  non- 
periodicity of  the  time  series,  and  two,  to  make  all  the  realizations 
the  same  length  of  336  hours  so  the  frequency  resolution  will  be 
equal  and  thereby  allow  ensemble  averaging. 

The  fourth  step  is  to  call  the  Fast  Fourier  Transform  (fT) 
algorithm.  The  basic  equation  in  the  algorithm  is 

V(n)  = Y v([j+l]At)e'^'’’j^^/'^'^^  (22) 

j=o 

where  V(n)  is  the  complex  Fourier  coefficient,  v([j+l]At)  is  the  value 
of  either  wind  speed  or  SO2  concentration  in  the  time  series,  and  NAt 
is  the  length  of  the  time  series.  A real  time  series  is  input  into 
the  FFT  subroutine  which  calculates  N/2  complex  Fourier  coefficients. 
The  (N/2+l)th  Fourier  coefficient  is  calculated  separately  and  in  this 


h[j] 


cos 


/ 


0 < j 
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analysis  contains  the  amplitude  at  the  frequency  0.5  cycles/hour. 
After  the  FFT  subroutine  has  provided  the  complex  Fourier 
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coefficients.  Equation  (23)  calculates  the  power  spectrum  with  an  ele- 
mentary band  width  of  0.00298  cycles/hour.  The  power  spectrum  is 


E^{n)  = ^|V(n)l2 


for  all  n from  1 to  N/2+1. 

Ihe  fifth  step  is  to  calculate  the  normalized  spectral  density 
by  dividing  the  spectrum  of  each  realization  by  the  variance  for  that 
realization. 

Finally,  each  frequency  is  ensemble  averaged.  The  result  is 
the  normalized  spectral  density  for  one  case  and  location. 
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Chanter  4 

I 

DATA 

i 

i 

Salt  Lake  Cit^,  Vicinity  | 

Of  the  twc  metropolitan  areas  selected  for  this  study,  the  | 

Salt  Lake  Valley  and  surrounding  areas  have  much  more  rugged  terrain,  | 

Figure  2.  Four  air  oollution  monitoring  sites  were  used.  The  sta-  | 

tions  are  located  at  Magna,  kearns.  Tooele  and  Bountiful.  1 

Sulfur  dioxide  (SO^)  is  the  effluent  monitored  in  this  study. 

It  was  chosen  because  of  its  long  half  life  and  the  relatively  high 
concentrations  at  the  four  sites.  The  major  single  SOo  source  in  the 
Salt  Lake  area  is  the  Kennecott  Copper  Corporation  smelter  located  at 
the  north  end  of  the  Oquirrh  Mountains  (kao  and  Taylor,  1977).  Since 
the  kennecott  smelter  is  a major  source  of  SOo.  especially  near  the 
Magna,  kearns  and  Tooele  sites,  it  is  assumed  to  be  a point  source. 

Also,  the  rate  of  effluent  ennssion  is  assumed  to  be  constant  through- 
out each  realization,  but  the  rate  may  vary  from  realization  to  real-  ; 

ization.  This  variance  is  of  no  consequence  since  deviations  from  the 
mean  or  trend  are  calculated. 

St.  Louis 

St.  Louis  is  a large  metropolitan  area  located  on  relatively 
flat  terrain  near  the  Mississippi  River  just  south  of  the  confluence 
of  the  Missouri  and  Mississippi  Rivers.  Figure  3 shows  the  monitoring 
stations  and  the  major  SC2  sources  affecting  the  stations. 
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Figure  2.  Base  map  of  the  Salt  Lake  Valley  and  Vicinity. 


Figure  3.  Base  map  of  St.  Louis. 


With  the  high  number  of  SO^  sources,  the  only  assumption  to 

C 

make  is  the  area  outlined  in  Figure  3 be  considered  an  area  source  and 
contains  the  majority  of  SO^  effluent  (Littman  et  al.,  1976).  Again, 
the  rate  of  effluent  emission  is  assumed  to  be  constant  throughout 
each  realization. 


Monitoring  Equipment 

The  Salt  Lake  City  vicinity  stations  are  operated  and  main- 
tained by  the  Utah  State  Division  of  Health,  Bureau  of  Environmental 
Health-Air  Quality  Section.  The  St.  Louis  stations  are  operated  and 
maintained  by  the  Environmental  Protection  Agency. 

The  anemometers  used  at  all  sites  in  the  Salt  Lake  City  vicin- 
ity and  St.  Louis  were  standard  rotating  cup  anemometers.  The  make  of 
the  instruments  varied,  but  all  have  a threshold  of  approximately  one 
meter  per  second  with  an  instantaneous  error  of  ±0.5  meters  per  second. 

The  SO^  monitors  vary  in  type  and  method  of  operation,  but  all 
have  a threshold  of  0.005  parts  SO2  per  million  (ppm)  and  an  accuracy 
of  ±0.005  ppm.  The  types  and  heights  above  the  surface,  of  the  anemom- 
eters and  SO2  monitors  follow; 

Magna;  Type  SO2  monitor:  flame  photometric. 

Height  above  surface:  30- foot  tower. 

Kearns:  Type  SO2  monitor:  conductimetric. 

Height  above  surface:  24-foot  tower  on  a one-story 
building. 

Tooele:  Type  SO2  monitor:  conductimetric. 

Height  above  surface:  24-foot  tower  on  a one-story 
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Bountiful:  Type  SO^  monitor:  conductrimetric. 

Height  above  surface:  24-foot  tower  on  a one-story 
building. 

St.  Louis:  Type  SO2  monitor:  all  flame  photometric. 

Height  above  surface:  30-meter  free-standing  tower  at 
all  sites  except  Site  8 which  has  a 10-meter  free- 
standing tower. 
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Averaging 

All  the  data  used  in  this  study  is  hourly  averaged.  Salt  Lake 
City  vicinity  station  averages  are  obtained  by  polling  the  instanta- 
neous wind  speed  or  SO2  concentration  once  every  six  minutes,  aver- 
aging the  values,  and  rounding  to  the  nearest  whole  mph  or  0.01  ppm, 
respectively.  The  St.  Louis  station  averages  are  obtained  by  polling 
the  instantaneous  wind  speed  or  SO2  concentration  once  every  minute, 
averaging  the  values,  and  rounding  to  the  nearest  0.01  m sec  or  0.0001 
ppm. 

Pasquill  (1962)  explains  the  effect  of  averaging  on  the 
spectra.  Two  problems  arise  in  the  averaging  process.  The  first  is 
to  reduce  the  fraction  of  the  total  variance  that  appears  in  the 
higher  frequencies.  The  second  is  to  introduce  aliasing. 

In  this  study,  the  fraction  of  total  variance  lost  in  the 
higher  frequencies  does  not  become  significant  except  for  frequencies 
greater  than  0.4  cycles/hour  (Table  1).  This  region  of  the  spectra  is 
small  compared  to  the  total  domain. 
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When  observations  are  at  discreet  intervals,  and  in  this  study 
the  interval  is  one  hour,  the  variations  associated  with  oscillations 


Table  1 


Effect  of  Hourly  Averaging  on  Spectral  Variance 


Frequency 

(cycles/hour) 

Fraction  of  Variance 
Remaining 

0.1 

0.97 

0.2 

0.88 

0.3 

0.80 

0.4 

0.62 

0.5 

0.45 

greater  than  0.5  cycles/hour  appear  at  a lower  frequency.  This  source 
of  error  is  called  aliasing.  If  aliasing  is  a problem,  the  spectra 
will  show  high  amplitudes  in  the  frequency  range  greater  than  0.3 
cycles/hour.  The  Bountiful  wind  speed  spectra,  but  no  other,  show  the 
high  amplitudes  and  the  reason  for  the  high  amplitudes  will  be  dis- 
cussed in  Chapter  6. 


Chapter  5 


STABILITY  CATEGORY  DETERMINATION 
Stability  Indices 

One  of  the  major  problems  in  the  air  pollution  meteorology 
field  is  determining  the  stability  of  the  atmosphere  in  the  surface 
and  planetary  boundary  layer.  Three  major  causes  of  this  problem  are 
diurnal  stability  fluctuations  at  a single  location,  horizontal  homo- 
geneity of  stability,  and  variations  in  stability  among  cases  that 
have  the  same  synoptic  situations.  Since  these  problems  exist, 
methods  of  classifying  periods  into  similar  groups  are  necessary  if 
any  analysis  of  SO2  and  wind  speed  fluctuations  are  to  yield  meaning- 
ful information. 

The  Salt  Lake  Valley  and  St.  Louis  have  different  topography 
and  thus  differing  stability  problems.  Because  the  valley  effect 
poses  unique  problems,  the  National  Weather  Service  Forecast  Office  in 
Salt  Lake  City  has  developed  a method  of  determining  stability  classi- 
fications for  the  Salt  Lake  Valley.  However,  St.  Louis  uses  nation- 
ally produced  air  quality  classification  products. 

The  differences  in  mixing  depths,  wind  speeds  and  air  pollu- 
tion potential  between  Salt  Lake  City  and  St.  Louis  was  found  by 
Holzworth  (1967).  He  found  that  throughout  the  year,  the  height  of 
the  mixing  depth  in  the  morning  averages  200  feet  lower  in  the  Salt 
Lake  Valley  than  at  St.  Louis.  The  average  wind  speed  in  the  morning 
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for  Salt  Lake  City  was  about  two  meters  per  second  slower  every  month 
than  at  St,  Louis.  The  afternoon  monthly  average  wind  speeds  were 
about  four  meters  per  second  slower  during  the  winter  months  at  Salt 
Lake  City  and  about  equal  during  the  summer.  The  air  pollution  poten^ 
tial  was  about  equal  on  the  average  even  though  the  Salt  Lake  City 
size  was  only  half  the  city  size  of  St.  Louis. 

Since  1967,  the  National  Weather  Service  Forecast  Office  at 
Salt  Lake  City  has  calculated  daily  an  Air  Stagnation  Index  (ASI) 
which  incorporates  the  height  of  the  mixing  depth,  mean  wind  speed 
through  the  mixing  depth  and  surface  temperature  (Jackman  and  Chapman, 
1977).  The  ASI  was  used  to  determine  the  stability  category  for  the 
Salt  Lake  and  adjacent  valleys. 

Seventy -eight  percent  of  the  stagnant  high  pressure  areas  in 
the  Salt  Lake  area  occur  during  the  months  of  November  through  Febru- 
ary (Jackman,  1968).  Since  the  stable  periods  are  of  most  interest, 
cases  were  only  selected  from  these  months. 

Jackman  and  Chapman  (1977)  found  conditions  in  the  Salt  Lake 
Valley  can  be  separated  into  three  categories  during  November  through 
February  (the  winter  months).  These  categories  are  ASI  < 200,  200  < 
ASI  < 1,000  and  ASI  > 1,000,  and  are  used  throughout  the  Salt  Lake 
Valley  analysis. 

The  first  category  ASI  < 200  is  called  the  stable  case.  This 
condition  occurs  when  cold  air  trapped  in  the  valleys  combined  with 
radiation  and  snow  cover,  results  in  strong  surface  inversions.  Usu- 
ally, warm  air  advection  above  the  inversion  tends  to  strengthen  the 
stable  condition  at  the  surface.  This  stagnant  layer  is  generally 
confined  to  below  6,000  feet.  Under  these  conditions,  diurnal 
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heating  is  unable  to  destroy  the  stable  layer  making  the  mixing  depth 

very  shallow.  Winds  below  the  inversion  are  very  light  and  often  | 

show  a diurnal  reversal  limiting  any  horizontal  transport  (Kao,  Lee 

and  Smidy,  1975).  I 

The  second  category,  200  < ASI  < 1,000  is  called  the  neutral 
case.  In  this  range,  synoptic  conditions  are  generally  characterized 
by  a subsidence  inversion  or  stable  layer  between  about  6,000  and 
12,000  feet.  Surface  heating  usually  allows  mixing  to  the  base  of  this 
stable  layer  which  gives  a moderate  mixing  depth  in  the  valley.  How- 
ever, the  base  of  the  stable  layer  may  be  at  or  just  above  higher 
mountain  areas,  and  may  severely  restrict  the  vertical  transport  of 
pollutants. 

The  third  category,  ASI  > 1,000  is  called  the  unstable  case. 

The  vertical  temperature  lapse  rate  approaches  the  dry  adiabatic  rate. 

Mixing  depths  can  extend  to  great  heights,  especially  in  the  presence 
of  approaching  frontal  systems. 

The  category  determination  for  St.  Louis  was  not  as  quantita- 
tive as  the  Salt  Lake  City  categories  because  no  localized  pollution 
index  was  available.  The  data  was  divided  into  two  categories,  stable 
and  neutral.  No  unstable  category  was  used  because  under  the  unstable 
conditions,  most  of  the  SO2  concentrations  were  zero.  Data  was  placed 
in  the  stable  category  if  St.  Louis  was  under  a stagnant  high  pressure 
area,  the  vast  majority  of  the  hourly  wind  speeds  were  below  three 
meters  per  second  and  the  hourly  wind  directions  were  between  090°  and 
180°.  Data  was  placed  in  the  neutral  category  if  no  stagnant  high 
pressure  was  over  the  St.  Louis  area,  no  precipitation  was  recorded, 
the  vast  majority  of  the  hourly  wind  speeds  were  between  three  and  six 
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meters  per  second,  and  the  hourly  wind  directions  were  between  090° 
and  180°. 


Quantity  of  Data 

I Data  for  the  Salt  Lake  Valley  cases  was  taken  from  1970  to 

I 1977  during  the  months  of  November  through  February.  If  a period  of 

[ three  consecutive  days  or  longer  with  the  same  stability  category  was 

i found,  then  the  period  qualified  as  a realization  for  the  respective 

• 

I stability  category. 

At  Magna,  both  SO^  and  wind  speed  data  was  used.  For  Magna, 
stable,  23  realizations  were  used  with  lengths  between  three  and  13 
I days.  For  the  neutral  case,  29  realizations  were  taken  ranging  from 

[ three  to  eight  days.  For  the  unstable  case,  17  realizations  were  found 

I varying  in  length  from  three  to  six  days. 

I At  Kearns  and  Tooele,  only  SO2  data  was  available  during 

I • stable  periods.  Eighteen  cases  at  Kearns  and  11  cases  at  Tooele  were 

f found.  All  cases  varied  in  length  from  three  to  13  days. 

I Bountiful  had  sufficient  data  to  calculate  both  stable  and 

1 neutral  wind  speed  and  SO^  spectral  densities.  For  the  stable  case, 

18  realizations  with  lengths  of  three  to  13  days,  and  for  the  neutral 
case,  10  realizations  of  length  three  to  five  days  were  found. 

The  data  used  in  the  St.  Louis  calculations  varies  from  the 
^ Salt  Lake  data  in  some  respects.  Since  only  1976  data  was  available, 

I realizations  were  taken  from  any  month  of  the  year.  When  24  consecu- 


tive hours  or  more  were  found  with  wind  direction,  speed  and  synoptic 
conditions  meeting  the  above  stated  criteria,  the  period  qualified  for 
use. 
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The  seven  St.  Louis  stations  had  the  same  realizations.  Sta- 
tions 1,  4,  5 were  grouped  together  and  called  urban  stations  because 
of  their  location.  Stations  8,  13,  20  and  21  were  grouped  and  called 
suburban  stations.  For  the  stable  case,  nine  realizations  of  SO2  and 
wind  speed  were  used.  For  the  neutral  case,  16  realizations  were 
found.  All  realizations  varied  in  length  from  one  to  three  days. 

Since  the  realization  length  is  much  shorter  for  St.  Louis 
cases  compared  to  Salt  Lake  Valley  cases,  the  lower  frequency  waves 
in  the  St.  Louis  spectra  will  be  suppressed.  Tl.is  fact  can  be  seen 
by  noting  the  spectral  densities  of  both  areas  in  Chapter  6. 


Degrees  of  Freedom 

Rayner  (1971)  gives  a method  of  estimating  the  degrees  of 
freedom  for  each  realization.  The  following  all  affect  the  degrees 
of  freedom:  the  number  of  points  tapered  at  each  end  of  the  realiza- 
tion, G;  the  length  of  the  realization,  D;  the  elementary  bands,  n/2; 
and  the  number  of  elementary  bands  in  non-overlapping  groups,  B.  An 
equation  to  estimate  the  degrees  of  freedom,  v,  is 
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B(D-G) 

n/2 


(24) 


where  G is  10  percent  of  D. 

Since  the  degrees  of  freedom  for  each  realization  in  a sta- 
bility case  are  summed,  the  D and  G in  Equation  (24)  are  the  sum  of 
all  the  realization  lengths  and  tapered  points,  respectively. 

The  degrees  of  freedom  vary  for  the  same  spectrum  because  the 


number  of  elementary  bands  in  non-overlapping  groups  increase  with 
increasing  frequency.  Only  one  elementary  band  is  in  the  groups  from 
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0.003  to  0.065  cycles/hour,  two  bands  in  the  groups  from  0.065  to 
0.119  cycles/hour,  five  bands  in  the  groups  from  0.119  to  0.298  ! 

cycles/hour  and  ten  bands  in  the  groups  from  0.298  to  0.500  cycles/ 
hour.  I 

The  degrees  of  freedom  for  the  shortest  and  longest  time 
series  were  calculated.  For  the  St.  Louis  urban  case,  which  is  the 
shortest  time  series,  the  degrees  of  freedom  are  5.2  for  the  0.003  to 
0.065  cycles/hour  range  increasing  to  52  for  the  0.298  to  0.500 
cycles/hour  range.  For  the  Magna  neutral  case,  which  is  the  longest 
time  series,  the  degrees  of  freedom  are  15.7  for  the  0.003  to  0.065 
cycles/hour  range  increasing  to  157  for  the  0.298  to  0.500  cycles/ 
hour  range. 


Chapter  6 


RESULTS  AND  CONCLUSIONS 
Auto-correlograms 

Autocorrelation  coefficients  were  calculated  only  for  Magna 
since  the  spectra  yields  more  information  than  auto-correlograms. 
Some  interesting  conclusions  can  be  drawn  from  the  auto-correlo- 
grams,  however.  The  autocorrelation  coefficients  of  SO2  and  wind 
speed  under  three  thermal  stability  conditions  were  investigated 
for  their  similarities  and  differences. 

The  auto-correlogram  under  stable  conditions  using  the  24 
hour  moving  average  shows  very  similar  coefficients  between  SO2 
and  wind  speed  up  to  a lag  of  13  hours,  (Figure  4).  The  co- 
efficients of  SOo  damps  out  for  longer  time  lags,  but  that  of  the 
wind  speed  continues  to  show  relatively  high  values  of  0.3  and 
0.2  at  24  and  48  hours  respectively.  This  correlation  which 
shows  the  diurnal  effects  of  the  wind  speed  in  the  Salt  Lake 
Valley  during  stable  conditions  could  be  a consequence  of  the 
effect  of  lake-valley  circulation. 

In  the  neutral  case  the  similarity  between  SO2  and  wind 
speed  auto-correlograms  (Figure  5)  is  less  than  in  the  stable 
case.  In  the  time  lag  of  the  first  five  hours  the  area  under  the 
curve  is  greater  for  SOo  than  wind  speed.  No  diurnal  correlation 
shows  in  the  wind  speed.  This  is  due  to  wind  speeds  being  less 
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affected  by  diurnal  circulation  in  neutral  periods. 

With  unstable  conditions  less  similarity  exists  between 
SO^  and  wind  speed  auto-correlograms  than  in  neutral  conditions 
in  time  lags  less  than  five  hours,  (Figure  6).  Again,  no  diur- 
nal variation  is  apparent  in  the  wind  speed  correlation  due  to 
the  intense  vertical  mixing. 

A direct  relationship  exists  between  the  integral  area  for 
time  lag  from  zero  to  four  hours,  and  the  stability  category 
(Figure  7)  for  SO2  and  an  inverse  relationship  exists  for  wind 
speed.  During  the  stable  periods  SO2  concentrations  are  generally 
better  correlated  than  those  during  neutral  and  unstable  periods. 

With  wind  speed,  higher  correlations  show  in  the  unstable  periods 
than  in  the  neutral  or  stable  periods. 

Two  similarities  exist  between  all  auto-correlograms. 

All  autocorrelation  coefficients  go  to  zero  between  time  lag  of 
four  to  five  hours  when  using  the  24  hour  moving  average,  and  the 
negative  correlations  are  very  similar  for  all  auto-correlograms, 
(Figure  7). 

Two  methods  of  calculating  the  autocorrelation  coefficients 
were  employed,  the  24  hour  moving  average  and  eliminating  the 
linear  trend.  Each  method  has  its  usefulness.  By  comparing 
Figure  7 with  Figure  8 it  can  be  seen  that  using  the  24  hour 
moving  average  to  calculate  the  autocorrelation  coefficients 
gives  the  best  similarity  in  the  auto-correlograms  of  the  hourly 
averaged  data  among  stability  categories  for  both  wind  speed  and 
SO2.  If  only  the  linear  trend  is  eliminated,  then  the  differences 
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correlograra  of  Magna  after  eliminating  the  linear  trend  with 
speed  (top)  and  SO,  (bottom). 
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in  the  integral  time  scales  among  the  stability  categories  will  be 
larger.  This  makes  for  easier  comparison  among  auto-correlograms 
for  different  stabilites. 

Cross -correloqrams 

Cross-correlograms  can  provide  useful  information  about 
correlations  between  different  stations.  The  cross -correlogram 
in  Figure  9 shows  correlations  of  SO2  between  Magna  and  Kearns. 

The  maximum  correlation  of  0.42  occurs  at  a lag  of  one  hour  when 
the  V£'(t)  is  calculated  from  the  Magna  data  and  the  v^'(t+T)  is 
from  the  Kearns  data.  This  result  shows  that  often  a change  in 
concentration  at  Magna  will  appear  at  Kearns  one  hour  later.  This 
occurrence  is  expected  since  Kearns  is  further  down  wind  from  the 
major  pollution  source.  A second  interesting  correlation  appears 
at  a lag  of  26  and  27  hours.  This  shows  a probable  diurnal  rela- 
tionship. It  will  be  shown  in  the  discussion  of  the  SO2  spectra  in 
stable  periods  that  a diurnal  variation  exists. 

Comparison  of  SO^  and  Wind  Speed  Spectrum 

' 4.  '' 

The  most  important  part  of  this  study  is  analyzing  the  sim- 
ilarities of  SO2  and  the  wind  speed  spectra.  It  is  the  objective 
of  this  study  that  similarity  between  the  spectra  can  be  established 
so  that  the  spectral  density  in  varying  spectral  bands  of  the  SO2 
spectrum  can  be  predicted  from  the  spectral  density  in  the  wind 
speed  spectrum. 

Since  the  Magna  site  provided  sufficient  data  to  analyze  the 
spectra  of  SO2  and  wind  speed  for  all  three  stability  categories. 


correlogram  of  Magna  vs.  Keams  under  stable  conditions 
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we  shall  first  analyze  the  spectra  at  Magna. 

The  SO2  and  wind  speed  spectra  for  stable  periods  were  simi- 
lar in  general  appearance,  but  the  wind  speed  spectrum  has  more 
pronounced  maxima,  (Figure  10).  In  both  spectra  a broad  maximum 
appears  in  the  frequencies  from  0.006  to  0.009  cycles/hours  (peri- 
ods of  168  to  112  hours).  These  low  frequency  waves  are  possible 
because  12  of  the  23  realizations  are  five  days  or  longer  thereby 
allowing  a wave  of  longer  than  112  hours  to  appear.  This  maximum 
is  more  pronounced  in  the  SO2  spectrum  than  in  the  wind  speed  spec- 
trum because  the  buildup  of  background  SO2  occurs  over  the  course 
of  at  least  three  days  during  a stable  period.  Another  broad  maxi- 
mum occurs  in  both  spectra  between  0.025  to  0.05  cycles/hour  (pe- 
riod of  40  to  20  hours).  This  maximum  encompasses  the  24  hour  diur- 
nal period.  The  amplitude  of  the  24  hour  period  is  especially 
strong  in  the  wind  speed  spectrum  and  is  probably  due  to  the  lake- 
valley  circulation  in  stable  periods.  Also,  two  other  peaks  corre- 
spond in  both  spectra;  the  0.083  and  0.125  cycles/hour  peaks  (12  hour 
and  8 hour  periods).  The  12  hour  peak  is  probably  due  to  a combined 
mountain- valley  and  lake  breeze  effect  and  will  be  discussed  in  more 
detail  later.  As  with  the  other  prominent  peaks  the  12  and  8 hour 
peaks  are  higher  in  the  wind  speed  spectrum. 

Unlike  the  0.006  to  0.009  cycles/hour  range  where  the  SO2 
spectrum  had  a more  pronounced  maximum,  all  peaks  with  frequencies 
at  higher  than  0.04167  cycles/hour  (24  hour  period)  show  more  pro- 
nounced maxima  in  the  wind  speed  spectrum.  This  fact  suggests  that 
a different  mechanism  is  affecting  waves  of  period  greater  than  24 
hours  compared  to  less  than  24  hours.  Since  the  wind  speed  spectrum 
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has  the  higher  peaks  in  the  frequencies  equal  to  or  greater  than 
0.04167  cycles/hour,  it  seems  likely  that  in  this  frequency  range 
the  wind  speed  spectrum  greatly  affects  the  spectrum  of  SOo-  But 
in  the  frequencies  less  than  0.04167  cycles/hour,  the  SO,  spectrum 

(to 

is  probably  influenced  more  by  the  builduo  of  SO2  trapped  in  the 
valley  than  by  the  wind  speed. 

The  neutral  and  unstable  spectra  can  now  add  more  substan- 
tiation to  the  above  claim,  (Figures  11  and  12).  For  the  neutral 
spectrum  the  wind  speed  shows  a 24  and  12  hour  maximum,  but  not  as 
predominant  as  the  wind  in  the  stable  category.  The  SO^  neutral 
spectrum  shows  a weak  maximum  around  24  hours,  but  no  12  hour  max- 
imum. Both  neutral  SO2  and  wind  speed  spectra  have  a broad  maximum 
centered  around  0.01  cycles/hour.  In  the  unstable  case  at  Magna 
a weak  peak  appears  around  24  hours,  but  no  other  significant 
peaks  are  evident.  For  SO2  small  peaks  appear  at  24  and  48  hours, 

P ’ but  at  no  other  frequency. 

In  the  stable  category,  the  wind  speed  and  SO,  spectra 
compare  similarily,  but  the  similarity  decreases  with  increasing 
atmospheric  instability  suggesting  that  variables  such  as  vertical 
motion  and  precipitation  have  more  effect  on  the  SO2  spectrum  than 
the  horizontal  wind  speed. 

This  decrease  in  similarity  between  the  SO,  and  wind  speed 
spectra  is  also  seen  in  the  St.  Louis  urban  spectra,  (Figures  13 
thru  16).  The  peak  in  the  urban  stable  wind  speed  spectrum  at 
0.095  cycles/hour  matches  the  SO2  spectrum  peak  at  0.083  cycles/ 
hour  (12  hour  period).  The  urban  neutral  wind  speed  has  no  peak  in 
the  0.08  to  0.1  cycles/hour  range,  but  a peak  at  0.083  cycles/hour 
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does  appear  in  the  urban  neutral  SO^  spectrum,  (Figure  ia). 


From  the  above  analysis  it  seem"  that  forecasting  the  SO2  spec- 
trum from  the  wind  spectrum  for  the  stable  category  can  be  successful 

I 

t ■ for  the  0.02  to  0.05  and  0.08  to  0.10  cycles/hour  ranges  (24  and  12 

hour  periods)  and  possibly  for  the  8 hour  period  and  0.003  to  0.015 
cycles/hour  ranges.  The  ability  to  forecast  the  SO2  spectra  for 
neutral  and  unstable  categories  is  doubtful. 

The  Inertial  Sub-range 

The  energy  decrease  in  the  spectrum  provides  information  about 
the  mechanisms  affecting  the  distribution  of  energy  under  various 
stability  conditions.  Pasquill  (1962)  stated  that  convectively  in- 
duced turbulence  in  the  inertial  sub-range  will  cause  energy  to  de- 
crease with  increasing  frequency  and  the  spectral  density  at  frequency 
n takes  the  form  of 

F(n)  « (25) 

where  e is  rate  of  energy  dissipation  per  unit  mass  of  fluid.  So  as 
n increases  in  the  inertial  sub-range  where  convectively  induced  tur- 
bulence occurs,  F(n)  will  decrease  at  a rate  proportional  to 
Thus,  Pasquills'  equation  is  named  the  "minus  5/3  power  law".  This 
means  the  decreasing  portion  of  each  spectrum  should  have  a slope  of 
-5/3.  Pasquill  also  stated  that  convectively  induced  turbulence 
extends  the  region  in  which  energy  decreases  rapidly  with  increasing 
frequency,  and  thereby  the  range  over  which  a "minus  5/3  power  law" 
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appl ies. 

To  evaluate  the  accuracy  of  this  theory  to  the  spectral 
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analysis  conducted  In  this  study,  a best  fit  straight! ine  was  esti-  i 

1 I 

mated  for  the  decreasing  portion  of  each  spectrum.  See  Tables  2 and  3 | I 

for  the  results.  From  Tables  2 and  3,  it  is  obvious  that  the  slope  I 


does  not  obey  the  minus  5/3  power  law  in  all  cases.  However,  the  in- 
ertial sub-range  is  extended  to  lower  frequencies  with  increasing  in- 


stability. This  is  shown  in  the  table  and  by  observing  the  respective  \ 

f 

spectra.  | 


Another  result  of  the  wind  speed  spectra  is  the  increase  in 
the  absolute  value  of  the  slope  with  decreasing  stability  for  stations 
in  the  Salt  Lake  area.  The  likely  reason  for  the  change  in  slope  is 
that  under  stable  conditions  the  large  scale  motions  are  suppressed  by 
the  inversion,  but  the  small  scale  motions  are  not  influenced  by  the 
inversion.  More  energy  can  remain  in  the  higher  frequencies,  and  have 
a smaller  decreasing  absolute  slope.  Under  neutral  or  unstable  condi- 
tions, the  microscale  turbulence  is  more  locally  isotropic  than  under 
stable  conditions,  therefore,  little  constraint  to  the  large  scale 
PKJtion  exists  and  the  lower  frequency  waves  absorb  more  energy.  Waves 
of  higher  frequencies  show  relatively  less  energy,  and  the  spectrum 
shows  a steeper  slope. 

Kao,  Paegle  and  Nonnington  (1974)  found  that  in  the  Salt  Lake 
Valley,  higher  frequencies  have  a more  pronounced  contribution  to  the 
spectral  densities  below  the  mountain  tops,  particularly  close  to  the 
surface.  Therefore,  the  slope  would  be  less  in  the  Salt  Lake  Valley 
than  under  equal  stability  conditions  in  a flat  area  such  as  St.  Louis. 
Since  a direct  comparison  between  Salt  Lake  stations  and  St.  Louis 
stations  cannot  he  made  because  of  different  methods  of  determining 
stability  category,  a fimi  conclusion  cannot  be  drawn  about  the 


I 

I. 


n- 


Table  2 

The  Slope  of  the  Wind  Speed  Spectra 


Station 

Stability 

Frequency  Range 
(cycles/hour) 

Slope 

Magna 

Stable 

0.14 

to  0.5 

-1.16 

Neutral 

0.13 

to  0.5 

-1.40 

Unstable 

0.045 

to  0.5 

-1.52 

Bountiful 

Stable 

0.18 

to  0.28 

-1.49 

Neutral 

0.175 

to  0.275 

-1.82 

St.  Louis-Urban 

Stable 

0.18 

to  0.5 

-2.25 

Suburban 

Stable 

0.15 

to  0.5 

-1.95 

Urban 

Neutral 

0.08 

to  0.5 

-1.63 

Suburban 


Neutral 


0.10  to  0.5 


-1.66 
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Table  3 

The  Slope  of  the  SO2  Spectra 


Station 

Stability 

Frequency  Range 
(cycles/hour) 

Slope 

Magna 

Stable 

0.04 

to  0.13 

-1.20 

0.13 

to  0.50 

-1 .66 

Neutral 

0.04 

to  0.10 

-0.92 

0.10 

to  0.50 

-1 .44 

Unstable 

0.04 

to  0.09 

-1.94 

0.09 

to  0.20 

-0.37 

0.20 

to  0.50 

-1.43 

Kearns 

Stable 

0.05 

to  0.15 

-0.75 

0.15 

to  0.50 

-1.74 

Tooele 

Stable 

0.05 

to  0.15 

-0.93 

0.15 

to  0.50 

-1.95 

Bountiful 

Stable 

0.175 

to  0.50 

-1.73 

Neutral 

0.125 

to  0.50 

-1.73 

St.  Louis-Urban 

Stable 

0.12 

to  0.50 

-1.62 

Suburban 

0.16 

to  0.50 

-1.71 

Urban 

Neutral 

0.10 

to  0.50 

-1.42 

Suburban 

0.13 

to  0.50 

-1.55 

wr 
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"mountain  effect"  that  Kao,  Paegle  and  Normington  (1974)  found.  But 
their  findings  might  explain  why  the  slope  at  Magna  during  stable  peri- 
ods was  -1.16  and  at  St.  Louis  -2.25  for  urban  and  -1.95  for  suburbs 
in  stable  periods.  The  neutral  spectra  for  the  same  locations  show  a 
slope  of  -1.40  and  at  St.  Louis  about  -1.65  for  urban  and  suburban 
stations. 

The  slope  in  the  inertial  sub-range  of  the  SO2  spectra  for  all 
sites  decreases  with  decreasing  stability.  Also,  the  decrease  is  bro- 
ken into  two  or  three  slopes  in  the  Magna  spectra.  The  difference  in 
the  SO2  and  wind  speed  decrease  can  be  partly  explained  because  SO^ 
concentrations  which  are  a scalar  quantity,  are  dependent  on  the  wind. 


Lake/River  Effect 

The  12  hour  period  in  the  wind  speed  in  the  Salt  Lake  Valley 
due  to  the  lake-land  circulation  was  confirmed  by  Kao,  Lee  and  Smidy 
(1975).  In  this  study,  the  12  hour  period  appears,  but  with  varying 
degrees  depending  on  stability.  From  Figure  17,  one  can  see  the  de- 
crease in  amplitude  of  the  12  hour  period  from  stable  through  unstable 
categories.  In  the  SO2  spectra,  a prominent  12  hour  peak  occurs  only 
in  the  stable  category  (Figure  18). 

A river  effect  is  found  in  the  data  obtained  at  St.  Louis,  but 
to  a much  smaller  degree  than  the  Great  Salt  Lake  effect  in  the  Salt 
Lake  Val ley  as  shown  inFigures  13,14,  15,  16, 17,  and  ip . For  stable  peri- 
ods, the  12  hour  period  set  up  in  the  St.  Louis  area  appears  between 
0.095  to  0.12  cycles/hour  in  the  wind  speed  spectra.  In  stable  periods 
the  SO2  spectrum  for  the  urban  stations  shows  a greater  amplitude  in 
the  frequencies  near  0.083  cycles/hour  (12  hour  periods)  than  the  sub- 
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urban  stations.  The  neutral  periods  for  both  urban  and  suburban  sta- 
tions show  little  energy  in  the  12  hour  period.  From  the  above  obser-  I 

vations,  we  can  conclude  that  the  river  effect  does  not  exist  far  from 
the  river  banks  and  only  appears  in  stable  conditions.  ) 

Effects  of  Local  Wind  Patterns  on  the  Spectra 

Of  the  four  Salt  Lake  vicinity  stations,  only  two,  Magna  and 
Kearns,  are  in  the  same  valley  (Figure  2).  The  local  wind  patterns  in 
the  different  valleys  have  specific  effects  on  the  spectra. 

The  SO^  spectra  for  stable  conditions  at  Magna  and  Kearns  are 
very  similar  as  shown  in  Figures  18  and  19,  which  is  expected  since 
both  sites  are  located  in  the  western  part  of  the  Salt  Lake  Valley. 

The  Tooele  SO2  spectrum  under  stable  conditions  as  shown  in 
Figure  20,  shows  one  marked  difference  compared  to  spectra  at  Magna 
and. Kearns.  The  difference  is  the  24  hour  peak.  Tooele's  24  hour 
peak  is  much  larger  than  Magna  or  Kearns  peak  indicating  a difference 
in  the  type  of  diffusion  in  the  Tooele  Valley  compared  to  the  Salt  Lake 
Valley.  The  Salt  Lake  Valley  is  closed  to  the  south  by  a mountain 
boundary  whereas  the  Tooele  Valley  is  open  to  the  north  and  south.  A 
tunneling  effect  is  set  up  and  enhances  the  lake-valley  circulation. 

The  spectra  at  Bountiful  as  shown  in  Figures  21  and  22,  dif- 
fers considerably  from  any  of  the  other  Salt  Lake  vicinity  stations  in 
two  ways:  first,  the  highest  peak  in  the  stable  wind  speed  spectrum 
is  12,  not  24  hours;  second,  much  more  energy  is  contained  in  the 
higher  frequencies  of  the  wind  speed  spectra  at  Bountiful  under  both 
the  stable  and  neutral  conditions. 

The  two  differences  are  attributed  to  the  lake  breeze  coupled 


cyeiM/  hour 

Figure  21.  Normalized  spectral  density  of  the  wind  speed 
at  Bountiful  under  stable  conditions  (top) 
and  neutral  (bottom). 
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with  canyon  winds.  Several  deep  canyons  are  located  in  the  Wasatch 

Mountain  Range  just  to  the  east-south-east  (ESE)  of  Bountiful  and  the  I 

Great  Salt  Lake  is  several  miles  to  the  west  (Figure  2).  Between 

the  hours  of  10:00  p.m.  and  9:00  a.m.,  the  prevailing  winds  are  ESE  I 

out  of  the  canyons  at  Bountiful  compared  to  the  prevailing  winds  at 
Magna  out  of  the  south-south-west  (Figures  23  and  24).  With  a lake/ 
canyon  wind  regime,  the  12  hour  oeak  is  explained. 

The  additional  energy  in  the  higher  frequencies  of  the  Boun- 
tiful wind  spectra  (Figure  21)  is  caused  by  an  increase  of  energy  in 
turbulent  eddies  due  to  the  canyon  winds.  Significant  energy  appears 
in  the  periods  ranging  from  3.5  to  2 hours.  At  this  frequency  range, 
aliasing  has  probably  occurred  and  indicates  considerable  energy  is  in 
•frequencies  higher  than  0.5  cycles/hour.  The  relatively  large  amount 
of  energy  in  the  higher  frequencies  indicates  that  even  under  stable 
conditions,  canyon  winds  produce  high  amounts  of  turbulence. 

The  canyon  winds  do  not  affect  the  SO2  spectra  in  the  higher 
frequencies,  but  the  lake-canyon  wind  shift  does  develop  a very  prom- 
inent 12  and  24  hour  peak. 
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CHAPTER  7 

SUMMARY 

Regarding  auto-  and  cross-correlograms  it  is  found  from  this 

study: 

The  similarity  decreases  between  wind  speed  and  SO2  auto- 
correlograms  with  decreasing  stability. 

2.  The  first  zero  autocorrelation  coefficient  occurs  near 
the  time  lag  of  four  to  five  hours  for  both  wind  speed  and  SO^. 

3.  The  24  hour  moving  average  gives  the  best  similarity  in 
the  auto-correlograms  of  the  hour  averaged  data  among  stability 
categories  for  both  wind  speed  and  SO2. 

4.  If  only  the  linear  trend  is  eliminated,  then  the  differ- 
ences in  the  integral  areas  among  stability  categories  will  be  larger. 

5.  For  stable  periods  the  time  lag  of  one  hour  has  the  high- 
est correlation  in  the  cross-correlogram  of  SO2  for  Magna  vs.  Kearns. 

The  following  results  were  found  from  spectral  analysis: 

1.  The  Magna  wind  speed  spectrum  greatly  affects  the  SO2 
spectrum  for  periods  equal  to  or  less  than  24  hours  under  stable 
conditions. 

2.  Under  stable  conditions  for  periods  longer  than  24  hours 
the  SO2  spectrum  is  probably  more  influenced  by  the  buildup  of  SO2 
trapped  in  the  Salt  Lake  Valley  than  by  the  wind  speed  spectrum. 

3.  The  similarity  between  the  spectra  of  wind  speed  and  SO2 
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decreases  with  decreasing  stability, 

4.  Forecasting  the  SO^  spectrum  from  the  wind  speed  spectrum 
for  stable  periods  is  possible  for  the  ranges  of  0.02  to  0.05  and 
0.08  to  0.10  cycles/hour. 

5.  For  both  wind  speed  and  SO2  spectra  the  range  over  which 
the  energy  decreases  is  extended  to  lower  frequencies  with  decreasing 
stabil i ty. 

6.  In  stable  periods  the  large  scale  motions  are  suppressed 
by  the  inversion,  but  the  small  scale  motions  are  free  to  absorb 
energy.  In  unstable  periods  little  constraint  on  the  large  scale 
motion  exists  and  the  lower  frequency  waves  absorb  more  energy. 

7.  The  slope  of  the  energy  decrease  at  Magna  for  the  wind  stieed 
spectrum  is  less  than  at  St.  Louis  in  stable  periods  probably  because 
in  mountainous  terrain  the  higher  frequencies  have  more  contribution 
to  the  total  energy  than  in  flat  terrain. 

8.  The  lake  effect  in  the  Salt  Lake  vicinity  produces  a 12 
hour  peak  in  the  wind  speed  spectra.  To  a lesser  degree  the  Missis- 
sippi River  produces  a river  effect  close  to  the  river  banks,  but 
only  in  stable  periods. 

9.  The  local  wind  patterns  in  the  different  Salt  Lake  vicin- 
ity valleys  have  specific  effects  on  the  spectra. 


APPENDIX  A 


The  first  program  calculates  the  autocorrelation  coefficients 
at  one  station,  and  the  second  program  calculates  the  crosscorrelation 
coefficients  between  two  stations.  Both  eliminate  the  linear  trend 


using  the  subroutine  Linear  Trend  Removal  used  in  the  first  program  of 
Appendix  B.  When  correlations  using  a 24  hour  moving  average  are  de- 
sired, replace  the  Linear  Trend  Removal  subroutine  with  Equation  (18) 
of  Chapter  3.  All  programs  in  Appendix  A and  B are  written  by  the 
author,  and  ICASE  is  the  number  of  realizations  for  each  case  and  M is 
the  length  of  the  respective  realization. 
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: PKOGRAK  TO  CAi.CULATE  THE  AUTOCORRELATION  COEFFICIENTS  FOR 
c ONE  station  anl  one  STASILITT  COIOITION 

OIMENSION  A(33b)  iXPRlNo.  (336)  tRTAU  (29.46)  • XPR1S2  (29>4»e>  . 

•TRTAU(46).  TXPRIS(46) 

ICASE  s 16 
KIS  1 

1 DC  9 KSl.ICASE 

IF  (N  .lE.  4)  «C  TO  101 
IF  (K  .lE.  6)  60  TO  1C2 
IF  (K  .lX.  10)  60  To  103 

IF  (K  ,lE.  13)  60  TO  104 

IF  (N  .Efi.  14)  60  To  105 

IF  (K  .EO.  15)  60  To  106 

IF  (N  .EC.  16)  60  TO  107 

1 01  96 

6C  TO  20C 

102  4S  12C 

60  TC  200 

103  k:  144 

6C  TC  200 

104  MS  166 

60  TC  200 

105  KS  192 

GO  TO  200 

106  MS  24C 

50  TO  200 

107  MS  336 
200  CONTINUE 

DO  10  is  1.M.12 

READ  (5.501)  A(I).  A(I'.l).  A(l42).  A(I43).  A(l44),  A(l45).  A(l46). 

• A(1-»7).  A(I46).  A(l49).  A(I410).  A(l4ll) 

10  CONTINUE 

501  FORMAT  (12F5.2) 

MA-.i2s  M - 12 
MM15S  w - 13 
CALL  LTR  (M.A.AV) 

00  13  LS  13.MM12 
13  XPKlMi(L)  S ACL) 

0 CALCIJ..ATE  THE  AUTOCORRELATION  COEFFICIENTS  USING  THE 
0 bwX  AND  UENKIKS  KETHOL 
00  15  is  1.46 
XPKISF-  0. 

XEOXJNS  0. 

RTAUCK.DS  C. 

XPR1S2(K.I)S  C. 

UC  16  wS  13.MM13 

AUTOCc  s XPRIM£(L)  • XPR1mE(L4:) 

XPRlS(ii  s XPRIME(L)4*2. 

XPRIL6  S XPRlME(L4l)*»2. 

RTAUCN.I)  s AUTOCO  4 RTAU(K.l) 

XPRISA  s XPRISG  4 XPRISM 
lb  XbOXJN  s XPR1L6  4 XSOAJK 

xpRisacK.i)  s (xpr:sm*xboxjk)«»o.5 

RXE  s KTAU(K.1)/XPRIS2(I..I) 

15  MM13S  Hm13  - 1 
9 COi;T1NU£ 

»r:TE  C6.62C) 

62C  FORMAT 1//10X, 'TOTau  AUTOCORRELATIONS  FOR  ONE  CLASS  OF  CLEARING 

• INDE.x'/') 
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)|  ■ 1 

H 00  20  Is  Ir«6 

TRTAUtJjs  0.  , 

TXPRISUJS  0.  I 

00  21  NS  1>!CASE 
TRTAUdl  S TRTAUd)  ♦ RTAUJK.I) 

21  TXPRIS(I)  s TXPRXStI)  ♦ XPRrS2(K.:J 

TOTRXE  s TRTA0(I)/TXPRIS(1)  \ 

RRITE  (6>62I)  I*  TOTRXE.  TRTAUd),  TXPRiSd) 

621  FORMAT  d0x.*TAUa',I3.5X.*TOTRXE*’.F10.7,5X»’TRTAU=*.El2.7.5X,* 

•TXPRISs. »EI2.7) 

20  continue 

KIs  KI  fl 

IF  «KI  ,E0.  2)  60  TO  I 

STOP  , 

SNl/  I 


H 
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C PR06RAM  TO  calculate  CROSSCORRELATION  COEFPECIENTS  BETWEEN  ANY 
C Two  STATIONS 

dimension  a(S36) iXPRlME <3S6>fRTAU (29,46) .XPRZSa(29,4a) • 

•TRT AU ( 46 ) , TXPR I S ( 46 ) • 6 ( 336 ) , B^R IME ( 336 ) . 

•CRT AU (29,46), TCTAU (46 ) , TCPRIS ( 46  > , CPRIS2 (29 ,46 ) • 

•R2ERC0(29),  RZERS6(29).  R2ERL6(29) 

XCASE  s 16 
DO  9 Xsl,lCASE 
ZF  (K  .LE.  A)  60  TO  101 
IF  (K  .LE.  6)  60  TO  102 
IF  (K  .LE.  10)  60  TO  103 

IF  (K  .LE.  13)  60  TO  lOA 

IF  (K  .£Q.  14)  60  TO  105 

IF  (K  .EG.  IS)  60  TO  106 

IF  (K  .EG.  16)  60  TC  107 

101  Ms  96 

60  TC  200 

102  Ms  120 

60  TO  200 

103  MS  144 

60  TC  200 

104  Ms  166 

60  TC  200 

105  Ms  192 

60  TO  200 

106  Ms  240 

60  TO  200 

107  Ms  336 

200  continue 

DO  10  IS  1,M,12 

READ  (5, SCI)  A(I).  A(l4l),  A(I«2),  Allwl),  A(It^),  AdwS),  A(l46), 
•A(l47),  A(I46),  AOO),  A(ITIO),  A(I411) 

10  continue 
501  format  (12F5.2) 

DO  50  IS  1.M,12 

READ  (5,501)  B(I),6(l4l), 6(142), 6(I«3),B(l44).b(l4S),B(l46), 

■B( 147), 6 (148) ,6(l49),B(I«10) ,8(1411) 

50  continue 

MMlZs  M • i2 
Mmi3s  M • i3 
CALL  LTR  (M,A,AV) 

CALL  LTR  (M,6,AV) 

DO  13  LS  13,MM12 
BPR!ME(L)  s B(L) 

13  XPRIMs(L)  s A(L) 

R2£RC0(K)s  0. 

R2£RSG(K)s  0. 

R2ERL6(K)s  0. 

DO  17  LS  13,MM12 

2ER0C0  s XPRXME(L)*8PRIME(L) 

2ER0SG  s XPRXME(L>«»2. 

ZEROLO  s BPRIME(L)«*2. 

RZERCO(K)  s ZEROCO  4 RZERCO(K) 

RZERSiii(X)  s ZEROSG  4 RZERSG(K) 

17  RZSRL6(X)  s ZEROL6  4 R2ERL6(K) 

DO  15  IS  1,46 
XPRISMs  0. 
xaox.^s  0. 

RTau(K,I)s  0. 

XPRIS2(X,I)S  0. 
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CPRIS^  = 0. 

CBOXJn  s 0. 

CRTAUlKd)  = 0. 

CPRZS2(K>I)  s 0. 

00  16  US 

AUTOCO  s XPRZMtiU)  • BPRIMECU-*!) 

CNCGCO  s BPRIM£(L)«XPR1ME(U«Z) 

XPRZSO  s XpRIMe(U)»*2. 

XPRZU6  s BPRIM£(U«'Z)«*2. 

CNbGSC  S BPRIME(U)»2. 

CNEGUv  s XPRZME(U-pZ)**2. 

XPRZSH  s XPRXSG  -»  aPRZSP 
XBOXJK  s XPRZLu  * XBOXJK 
CPRZSix  s CnEGSO  ■»  CPRZ5H 
CBOXJK  s CnEGUG  « C60XJK 
RTAu(a«Z)  : AUTOCO  ♦ RTAU(K>Z) 
lb  crtau(kiZ)  s ckegco  * crtau(k>:) 

XPRZS2(KtI)  s (XPRZSM>XBOXJK)««0.5 
CPRZS2(K.Z)  s (CRRZSM*CBOXJK)«.0.5 
15  MMlSs  MM13  > 1 
9 CONTINUE 
T2ERCC  s 0, 

T2ERSG  s 0. 

TZERU6  s 0, 

00  22  Ks  I.ICASE 
TZERCO  : RZERCO(K)  <»  TZERCO 
rzCRSO  s RZERSG(K)  4^  TZERSO 
22  TZERLG  S RZERUG(K)  4.  TZERUG 
TOTZCl  s TZERCO/TZeRSO 
T0TZC2  s TZERCO/TZERUG 
TOTZC3  5 TZERCO/ (TZERSO-TZERLG) ••0.5 
•kite  (6,622)  TOTZCl.  T0TZC2,  T0TZC3 
622  PORNAT  (lOx, 'TOTZCls* ,F7,b,5X, »T0TZC2s* ,P7.A,5X, 'TCTZCSr • ,P?,U//) 
mRZTE  (6,620) 

620  PORMAT(//10X, 'TOTAu  AUTOCORRELATIONS  FOR  ONE  CLASS  OF  CLEARING 
• ZNOEa'/) 

00  20  zs  i.ue 

TRTAU(Z)s  0, 

txpr:s(I)s  0. 

TcTAU(Z)  s 3. 

TCPRZS(Z)  s 0, 

00  21  l.ZCASE 

TRTAU(i)  s TRTAU(Z)  ♦ RTAU(K,I) 

TXPRIS(Z)  s TXPRIS(Z)  ♦ XPR1S2(K,Z) 

TCTAU(l)  s TCTAU(Z)  ♦ CRTAL(K,I> 

21  TCPRZS(Z)  s TCPRZSd)  • CPRIS2(K,Z) 

TOTHXE  s TRTaU(Z)/TXPRXS(I) 

TOTCXE  s TcTAU(I)/TcPR:S(1) 

WRZTE  (6,621)  I,  TOTRXE,  TRTaU(I),  TXPRIS(I),  TCTCXE, 

•TCTaU(Z),  TCPRZS(l) 

621  FORMAT  ( 'TaUs* . I3,3X, 'TOTRXEs* .FT.U.SX, ‘YRTAUs* ,E9. A, IX. 'TXPRISs* , 
•E9.A.1X. ’ToTCXEs* ,F7,a,3X, 'TCTaUs* ,E9.4.ix, 'TCPRlSs* .E9.a) 

20  CONTINUE 
STOP 
ENO 


APPENDIX  B 


' 
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SPECTRAL  DENSITY  PROGRAMS 

The  spectral  densities  are  calculated  for  the  Salt  Lake  Valley 
and  vicinity  stations  in  the  first  program  and  for  the  St.  Louis  urban 
and  suburban  stations  in  the  second  program.  The  two  subroutines. 
Linear  Trend  Removal  and  Cosine  Bell  Window,  are  modified  in  the  first 
program  so  that  the  12  first  and  last  data  points  are  not  used.  This 
step  is  necessary  since  12  extra  data  points  were  added  to  the  Salt 
Lake  Valley  and  vicinity  data  so  that  in  the  programs  of  Appendix  A 
the  ii^  and  the  u^^  of  the  24  hour  moving  average  could  be  calculated. 


j 
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C PROORAM  FOR  CALCULAT1N6  THE  POHER  SPECTRUM  PROM  A SINGLE  MONITOR 
c POR  MANY  pollution  CASES  USING  THE  DIRECT  TRanSPORM  METHOD  ON  PAGE 

C S2  • G2  RaTNER  (1971)  and  MODIFIED  SUBROUTINES  PROM  LEE  AND  HATCH. 

0 1 MENS I ON  A ( 336 ) • C ( 16B ) • I PK ( 1 049 ) » E ( 336  > . XPR IME ( 336 ) . 

•C V AR ( 29 ) . CSPEC (29.I69).SPEC(29.169)>XVAR(29). TSPEC ( 169 ) > 
•TCSPEC(169)f  SPAV6(169) 
complex  c«  «amn 
equivalence  (EiC) 

NS  336 
N2  s N/2 
N3  S N2  « I 
PS  O.I 
ICASEs  IS 
KIs  1 

I 00  9 Ksl.ICASE 

IP  (K  .LE.  6)  GO  TO  lOX 
IP  (K  .LE.  10)  GO  TO  102 

IP  (K  .LE.  11)  60  TO  103 

IP  (K  .LE.  lA)  GO  TO  106 

IP  (K  .LE.  IT)  GO  TO  110 

IP  (K  .EG.  16)  60  TO  107 

101  MS  96 

GO  TO  200 

102  MS  120 

GO  TO  200 

103  MS  ltt« 

60  TO  200 

106  MS  216 

GO  TO  200 
110  M s 290 

GO  TO  200 

107  MS  336 
200  CONTINUE 

12  s (n«M)/2  * 12 

MM12  s M • 12 

MN2A  s M • 24 

M24I21  s M • 24  •«-  12  1 

DO  10  IS  1.M.12 

READ  (5*501)  A(I)*  A(1<»1)*  • A(I4>3)»  A(I<44)>  A(l4>S)*  A(lf6). 

•Ad^T)*  A(14^G)*  A(I«9)*  A(l4-10).  Ado’ll) 

501  format  (12F5.2) 

10  CONTINUE 

CALL  LTR(M«A*AVG) 

CALL  RINOOk  (M*A*P) 

C GIVE  all  THE  periods  THE  SAME  LENGTH  ST  ADDING  EQUAL  NUMBER  OP 
C 2ER0S  TO  both  ENOS  OP  THE  DATA. 

DO  S3  IS  1,12 

53  XPRlME(I)s  0. 

DO  52  IS  1,MM24 

52  XPRIME(I4I2)  s A(l4l2) 

DO  54  is  M24I21,N 

54  XPRIME(I)s  0. 

c calculate  The  spectrum  using  the  past  pourier  transform  algorithmn. 

DO  45  KKs  1,N 
45  E(KK)  S XPRIME(KK) 

call  PPTR(E*GAMN>N,IMK) 

00  30  KKs  1,N2 
30  c(KK)  s c (kk) /Float (N2) 

GANNS  6AMN/PL0AT(N£) 


ri 
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DC  35  IS  I.N2 

35  SPEC(K.I)  s CABS(C(I) )*»2./2. 

SPEC(KiN3)  s CABS(&AMN)«*2. 

CVAR(K)  s 0. 

00  so  IS  Itl42 

50  CVAR(K)  s SPEC(KiI)  -p  CVAR(K) 

CVAR(R)  S CVAR(K)  -p  SPEC(K»N3) 

DO  21  IS  l.NS 

21  CSPEC(K«1)  s SPeC(K«I)/CVAR(K) 
ys  0. 

DO  70  IS  1>N 
70  rs  Y •*  XPRIME(I)**2. 

XVAR(K)  S Y/N 
9 CONTINUE 

C ENSEMBLE  AVERASE  ThE  CASES  POR  EACH  STABILITY  CLASS. 

TCVAR  s 0. 

TXVAR  5 0, 

DO  22  KS  ItICASE 
TCVAR  s CVAR(K)  * TCVAR 

22  TXVAR  s XVAR(K)  TXVAR 

DC  23  IS  1.N3 

TCSPECd)  s 0. 

TSPECd)  s 0. 

DC  24  KS  1. I CASE 
TSPEC(l)  s SPEC(K,I)  ♦ TSPEC(l) 

24  TCSPECd)  s CSPEC(K.1>  -p  TCSPECd) 

23  CONTINUE 
DO  25  is  1>N3 

TSPECd)  s TSPECdJ/dCASE«TCVAR) 

25  TCSPECd-)  s TCSPECd >/ICASE 
C AVERAGE  The  POWER  SPECTRA  IN  THE  HIGHER  FREOUENClESt (0.065 
C CYCLES/HOUR  and  HIGHER). 

00  27  IS  22>40i2 

27  SPAV6(1*1)  s (TCSPECd)  -p  TCSPEC d-Pl) ) /2. 

00  26  IS  42.97.5  I 

26  SPAV6(l-p4)  s (TCSPECd)  -p  TCSPECd+l)  •*  TCSPECd-p2)  ♦ TCSPECd-pS) 

• ♦ TCSPEC(I-p4)  )/5. 

DO  29  Is  102.152. 10 

29  SPAV6(I^9)  S (TCSPECd)  ♦ TcSPECd  + l)  ♦ TCSPECd'P2)  ♦ 

• TCSPEcd*3)  ♦ TCSPEC(1^4)  ♦ TCSPECd+S)  ♦ TCSPECd-pG)  -p 

• TCSPECd-pT)  "p  TCSPEC(I-p6)  -p  TcSPEC(1-p9)  )/10. 

DO  31  IS  162.166 

31  SPAV6(166)  s SPAV6(166)  -p  TCSPECd) 

SPAV6(166)  s SPaV6(166)/7. 

■RXTE(6.610) 

610  format  {//lOX. 'WAVE’ .6X. *FREQ' .8X. 'HOURS* .9X. 'DAYS' .7X. 'SPECTRUM* , 

•9X .' C-SPEC ’. 9X .' SLOPE '/ ) 

IMOs  0 

WRITE  (6.612)  IMO.  TSPEC(l).  TCSPECd) 
format  (10X.I4.62X.E10. 5. SX.ElO. 5) 

DO  26  is  2.N3 
IMXs  I • 1 

hrss  Float (N) /float (iMi) 

DATS  S HRS/24. 

FREO  s 1,/hRS 

WRITE  (6.611)  IMI. FREG. HRS. DAYS. TSPECd). TCSPECd). SPAVGd) 
format  (10X.I4.5X.F6.S.SX.F9.2.SX.F6.2.5X.E10.5.5X.E10.5.5X.E10.5) 

CONTINUE 

WRITE  (6.621)TXVAR.  TCVAR 

format  {/lOX. 'variance* .ex. E1Q.5.5X. 'SPEC  VARIANCE' »2X»E10. 5////) 

Kl  S RI  ♦ 1 

IF  (KI  .EG.  2)  60  TO  1 
STOP 
END 


612 


611 

26 

621 
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Cm*>  tlNEAR  TKEnO  removal  RROORAM 
SUBROUTINE  LTR(N»X«SX) 

DIMENSION  X(N) 
fNsN-<« 

MMAteN«12 

SXsO.  \ 

SM*0.  ] 

00  II  instMMU  j 

II  sxssx^xd)  j 

SXSSX/FN 
00  a 1SI0,MM12 
X(2)SX(1>.SX 
13  SPsSP^XCD-FLOATd} 

AVIs(FN*i.)/E. 

AN6s-I2.*SP/(FUb(FN«<»2-I)  ] 

00  IS  ISI3>MMI2 
FIsI 

15  X(I)aX(I)-AN6*(AVI-Fl) 

RETURN 

ENO 


C»«*»  COSINE  BELL  FINDOf 

SUBROUTINE  WlNOOH(N«XtP) 
dimension  X{N) 

«s:PbNbi2, 

NTENTmxG 

PlftsS. 14159265/ (0-12. ) 
NPSN*1 

00  2 IsI3»NTENTh 

FIsI-12 

jNPaNP-l 

C0SaELs0.S«(I.-C0S(Fl*PI6) ) 
X(I)bX(I)«COSBEL 
2 X(JNP)&X(JNP)«C0SBEL 
RETURN 
ENO 
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PROGRAM  FOR  CAUCULAUNG  SPE.CTRAu  DENSITITES  FR0»'  SEVERAL  MONITORS 

pOR  many  pollution  cases  using  the  direct  transform  method  On  Page 
S2  - 63  RAYnER  (1971)  and  SUBROUTINES  FROM  uEE  AND  HATCH. 

DIMENSION  a(336)«C(166) •IWK(1049) »E(336) .XPRIM£(336) . 

•CVAR (29> « cSPEC (29*169) fSPEC (29,169) 'XVaR (29) , TSPEC (169) , 
•AVG(lC),TCSREC(10,le9).  SPaV6(10,169) , AVSPEC (169) ,TSPaVG(169) 
COMPLEX  C,  GAMl. 

EGUIVALENCE  (E,C) 

Ns  336 
N2  s N/2 
N3  s N2  'Y  1 
Ps  0.1 
ICASE:  9 
Kls  1 

1 DO  9 NslilcASE 


IF 

(K 

.£(t. 

1) 

Ms36 

IF 

(K 

.EQ. 

2) 

ms2a 

IF 

(K 

.EQ. 

3) 

HsAe 

IF 

(K 

.EG. 

A) 

MsAo 

IF 

(K 

.EG. 

S) 

Ms  6u 

IF 

(K 

.EG. 

6) 

Ms  36 

IF 

(K 

.EQ. 

7) 

Ms  2A 

IF 

(K 

.EG. 

6) 

Ms  2A 

IF 

(K 

.eg. 

9) 

HS  2A 

12 

s 

{N-M)/2 

LaSTZs  N-I2+1 

DO 

10 

IS  1 

12 

READ  (S,501)  A(I),  Ad^l),  A(I.,2),  A(I-y3)>  Ad-rU).  Ad<*6).  Ad«6), 
•Ad<»7},  A(lY-e),  Ad4-9),  A(l<»10),  Ad-Yll) 

SOI  format  (12FS.2) 

10  continue  ' 

CALL  LTR(M,A,AVG) 

Call  <i1NOOi»  (m,a,P) 

C GIVE  ALL  the  periods  THE  SAME  LENGTH  BY  ADDING  EQUAL  NUMBER  OF 
C ZEROS  TO  BOTH  ENOS  OF  THE  DATA. 


DO 

53  IS 

1.12 

S3 

XPRIMEd) 

s 0. 

DO 

52  IS 

1,M 

52 

XPRIMEd* 

12)  s 

Ob 

5A  Is 

LaSTZ 

5A 

XPRIMEd) 

s 0. 

c calculate  The  spectrum  using  the  fast  fourier  transform  algorithmn. 

DO  AS  KKs  1,N 
AS  E(KK)  s XPRIME(KX) 

CALL  FFTR(E,6AMN,N,1WK) 

DC  30  KKs  1,N2 
30  C(KK)  s c ( KK] /Float (N2) 

GAMNS  GAMN/FL0AT(N2) 

DO  3S  IS  1,N2 

3S  SPEC(K,X)  s CABS(C(1))**2./2. 

SPEC(K,N3)  S CABS(GAMN>a«2. 

CVAR(K)  s 0. 

DO  so  IS  1,N2 

50  CVAR(K)  s SPEC(K,1)  a CVAR(K) 

CVAR(K)  s CVAR(K)  ♦ SPEC(K,N3) 

DC  21  is  1,N3 

21  CSPCC(K,I)  s SPEC(K.I)/CVAR(K) 

Ts  0. 

DO  TO  Is  l.N 
70  Ts  Y ••  XPRlMEd)»»2. 

XVAR(K)  S Y/N 
9 CONTINUE 


1 
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E.<sEMaL£  averase  The  ;ases  «fOR  each  stabilitit  class. 

TCVAR  S C. 

TXvAR  s 0. 

OC  22  R:  A.ICASE 
TCVAR  s CVaRiK)  ♦ TCVAR 

22  TXVAR  s XVaR(KI  * TXVaR 
DO  22  Js  1,N2 
TSPECCI)  s 0. 

DO  2R  K=  l.ICASE 
TS?EC«I»  s SPEC(K,I) 

24  TcSPEClKI.;)  s CSPECCK.I) 

23  CONTINUE 


♦ TSHECd) 

tcspec(k:i 


1) 


DO  25  is  1,N3 

TSPECd)  s TSPfcC(i)/(ICAS£»TCV,-.R) 

25  TcS'>EC(KI,l)  s TCSPEC  (Kl , 1 ) /ICaSE 

; AVERauE  TnE  pokier  spectra  in  the  hISHEB  PREOUENCIES.  10.065 
C CTCLES/HOUR  AnC  H16HER). 

DC  27  1=  22.40.2 

27  SPAV6  IK:  . 1*1  »s  (TCSPEC  (Kl , 1 )*TCSPEC  (K:  . l-^l ) ) /2. 

DC  26  Is  42.97,5 

26  SPAV6lNl,l*4)s(TCSPEC(Kl,l)*TCSPEC(KI,l-.i)*TrSP£CIM.l42)4 
• TCSPEC(KI,1-»3)*TCSPEC(KI.I4'4)  1/5. 

DO  29  is  1o2.152>1C 

29  SPAV6{k:,I^9)s{TCSPECCKI.I»‘*TCSPEC(K1,1^1)*7CSPEC«KI.1'*'2)4 
-TcSPEC(KI,i*3)*TCSPEC<Kl.l*4)*7CSPEC(KI,l45)4TCSPEC(KI,l4.6)-t. 
•TcSPEC  IKI , !♦?  )+TCSPEC  (Kl , 14.6  )*TCSPEC  (K1 . 14.9)  J /lO . 

00  31  is  162.166 

31  SPAV6IK!.l6e)sSPAV0(K1.16e)4TCSPEC(K1.2) 
SPAV6(Kl.l6e)sSPAV0(K1.16B)/7. 

TaV6s  0. 

DO  el  USl.lCASE 
61  TAVSs  TAV6..AV6U) 

•RITE  (6,632)  TAV6 
632  PORMaT  (10x»P10,5///> 

R1  s M ♦ 1 

IF  (Kl  .LE.  7)  60  TO  1 

C AVERAGE  EACH  PREOUENCT  OF  THE  THREE  URBAN  STATIONS 
DC  62  is  2,N3 
DO  63  US  1,3 

TSPAV6(I)s  TSPAVGCI)  .»  SPAVG(J,n 
63  AvS?EC(I)s  AVSPECd)  ♦ TCSPEC(w.I) 

TSPAV6(1)5  TSPAVGCI )/3, 

AVSPECd )s  AVSPECd )/3. 

iMlS  I-i 


PRESS  FLOAT (lHl)/FtaOAT(N) 

62  PRITE  (6.633)  IMI.  PREO>  AVSPECd)*  TSPAVGd) 

c average  Each  freouenct  of  the  four  suburban  stations 

DO  72  is  2,N3 
DO  73  js  4,7 

TSPAVGcDs  TSPAVGCI)  * SPAVGCJ.X) 

73  AVSPECCDs  AVSPECd)  * TCSPEC (J. 1 1 
TSPaVGCDs  TSPAV6C1)/4. 

AVSPECCDs  AVSPECd )/4. 
iMls  1-1 

FREGs  FLOAT (iMi) /Float (N) 

72  WRITE  (6,633)  IMI,  FREOi  AVSPECCl),  TSPAVGCI) 

633  format  C10x,I4,F12.S,E15.5.E15.5) 

STOP 

END 


1 


C*««a  COSINE  sell  WJNqOm 

subroutine  »lNUO»r(i>.»X«P) 

DIMENSION  X(N) 

OcPaN 

NTENThsO 

P16sB.1«I5«265/G 

NRsNal 

00  a XbI»NTENTh 
Pill 

UNPsNP»I 

COSBELsO.Sad.-COSlFIaPIO)  I 
XdlsXtliaCOSBEL 
a X(JNP>BX(UNP)*COSb£L 
return 

ENC 


I 


I 


f. 
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